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1. Introduction 

The transverse Ising model (TIM) is described by a two-state Ising Hamiltonian with a 
term representing a field transverse to the spins, 

H = -J2J^JS^S]-QY.S!, (1) 

ij i 

where Jij > 0, Q is the transverse field, S'f and S"? are Pauli spin-| operators and 
the first sum is over the nearest neighbors spins on the lattice. The model has been 
firstly applied to describe the phase transitions and the properties of hydrogen bonded 
ferroelectrics [H |2] and magnetic ordered materials [3] . This model in one dimension has 
no phase transition at finite temperatures; however, at zero temperature it is ordered 
up to the critical value of the transverse field. The model has been solved exactly in one 
dimension [H El |6]. In high dimensions there are approximations for low-temperatures 
or high-temperatures regions [3 [8] . All other calculations are based on mean field type 
approximations. An effective field theory has been presented which improve over mean 
field results |9j. Since then many results have been obtained based on the effective 
field theory. More recently the model has been used to study the phase diagrams of 
nanowires systems [10] and magnetization of nanoparticles [11] . The objective of this 
paper is to present rigorous upper bounds for the critical couplings. We will apply the 
results for the d=2 square lattice and the d=3 cubic lattice. 



2. Generalized Callen's identity for the transverse Ising model 



In this section, we describe the methodology used by Sa Barreto et al [9j to derive 
an identity for the two-spin correlation function of the transverse Ising model. The 
procedure used in this deduction was presented in reference [9] to obtain an exact 
relation for the order parameter <Sf> which generalizes Callen's identity [12]. The 
longitudinal two-spin correlation function <SfSl> can be calculated from 

Tr{e-^^StSl) 



(2) 



Tr(e-/3^^) 

where H is given by ([1]). The Hamiltonian can be separated into two parts, H = Hi + H', 
where Hi includes all parts associated with site i and H' represents the rest of the 
Hamiltonian. A direct calculation leads to 



< sist >- 



Sf 



Tr^i^Sfe-^ 
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(3) 



where Tr(j) represents the partial trace with respect to site i and A = 1 — 
^-pH.^-pH' ^I3{H,+H') ^ Equation ^ is an exact relation. However, it is difficult to be 
used. Therefore, we will make an approximation based on the following decoupling, 



Sf 
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Inserting (jl]) into Q and using the fact that < A >< 1, we obtain, 

< StSt >< (St '''- ':' ■ (5) 



By expanding A we see that the approximation is correct to the order of Moreover, 
it is consistent with the apphcation of the correlation inequahties that will be used later 
to obtain the upper bounds for the critical couplings. In the next steps we will keep 
only the = sign of . 

Let us write Hi = —EiSl — ^S^, where Ei = J2j JijSj- Diagonalizing and taking 
the partial trace over i , we get for the longitudinal spin correlation function, < crfa^ >, 
where cxj = 2 Si, 

J2i Jij^j 



z\2 



tanh(/3^(2^])2 + (^J,,.a|)2)). (6) 

Introducing the exponential operator e'^°'^\f{x) = f{x + a) , D = we obtain, 
< afa^ >= (afe^^^^^^"^'^^) ■ f{x) 



\x=0 



3 

3 



where f(x) is given by 



X 



fix) = , tanh(/3) J(2fi)2 + x^ (8) 

'{2Qy + x^ ^ 



Note that f{x) = -f{-x). 

Expanding the exponential in ([7]) and considering that {crlY = 1, we obtain, 

< alat >= {af [] [cosh( J,,D) + sinh( J,,D)] ) ■ f{x) U=o (9) 

3 

By a similar procedure the transverse two-spin correlation function <afai> is obtained, 

< af af >= (af [] [ cosh( + sinh( J^jD)] ) ■ g{x) \^=o (10) 

3 

where g(x) is 

2fi 



g{x) = . tanh(/3)^(2fi)2 + x^ = g{-x) (11) 



The expectation value of cxf is given by, 

< >= (n [cosh(J,,D) + a,^sinh(J,,D)]) ■ (7(x) U=o (12) 
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3. Application to d=2 square lattice and d=3 cubic lattice. 



3.1. d=2 lattice 

Considering the four neighbours of i in ([9]), expanding the product, applying the 
exponential operators appearing in the powers of cosh(JjjD) and sinh(JijD) in f(x), 
we obtain, 

j j<k<m 

where 

A2 = ^[/(4J) + 2/(2J)]>0 
S2 = ^[/(4J)-2/(2J)]<0 

(14) 

and j, k and m are neighbours of i and /(..) is given by ([8]). 



3.2. d=3 lattice. 

After a similar calculation we obtain for the cubic lattice, 

j j<k<m 

+ Cs J2 <(^!^]^l«^p > (15) 

j <k<m<n<p 

where, 

^3 = ^[/(6J) + 4/(4J) + 5/(2J)] > 

S3 = ^[/(6J)-3/(2J)]<0 

^3 = ^[/(6J)-4/(4J) + 5/(2J)]>0 

(16) 

and j, k, m, n and p are neighbours of i and /(..) is given by ([8]). 



4. Exponential decay of the two-point functions and the upper bounds. 

Upper bounds for the critical temperature T^. for Ising and multi-component spin systems 
have been obtained by showing (for T>Tc) the exponential decay of the two-point 
function [131 El US]- The procedure to obtain these upper bounds for the critical 
couplings of the tranverse Ising model is the following: we start from a two-point 
correlation function equations f|T3l) and f|T5l) and we make use of Griffiths inequalities 
(Griffiths 1, 11) ITS! [TU] and Newman and Lebowitz inequalities pTf [TU] . A proof 
of Griffiths inequalities has been given for the XY model with no external field \16\ . 
Extensions of Griffiths-Kelly- Sherman inequalities to quantal systems, under external 
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fields, both longitudinal and transverse, have been proved [ISl[I2|. The physical reason 
why the Griffiths and similar inequalities are valid for the quantal XY-type Hamiltonian 
is that the off-diagonal interaction, namely, Hi{x) = J2a Ja^% {J a — 0); produce the 
decrease of the ferromagnetic correlation among the (T|-spins, but it is not sufficient 
big to create a cooperative effect to induce an antiferromagnetic correlation. In other 
words, one can say that Hi{x) is a dynamical random force acting on z-z correlations 
[Is] . We establish the inequality for the two-point function <afa^>, 



< olal >< J2 < >, < a, < 1, (17) 

i 

which when iterated [H] implies exponential decay for T > Tc. 
4-1- Upper bounds for d= 2. 

From Eq. (fT3l) . using Griffiths II (<(Tf (T|cr^cr^>><(Tf (T|><cr|(T.^>) in the second term 
and considering B2<0, we get, 

< (Tfcrt >< E «i < ^1 >' (18) 

j 

where j ^ k are neighbours of and 

Oj = A2-\B2\< a]al >id (19) 

4-2. Upper bounds for d=3. 

From Eq. (fT5|) . using Griffiths II in the second term -B3<0, Newman's inequality 
{<alF>< J2j <o'iaj><dF/daj>, F are polynomials with positive coefficients) 
combined with Griffiths I {<o'%>< 1) on the third term C3>0 , we get, 

<afat><Y.a,<afa]>, (20) 
j 

where j ^ k are neighbours of i, and 

Uj = As- \Bs\ < a]al >i, +5C3 (21) 

4-3. Numerical Results 

The two-spin correlation functions appearing in (fT9|) and (12T|) < c^a^ >id is the one- 
dimensional model two-spin correlation function separated by a distance of two lattices 
sites. For the one-dimensional transverse Ising model the exact value of this function at 
the critical value Qc = J is 

1(2.2 H\2) 

where H{n) = l"^^2'"^^...(n — 1). For the one-dimensional Ising model the exact value 
of the spin correlation function separated by a distance of two lattices sites is: 

< >id=tanh^/3J (23) 



<a]al>^d=-{-) 8^ (22) 
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In the following, we will use result in (HM and ([21]) to calculate the upper 
bound for at Tc = and result fl23l) in f|T9|) and fl2Tl) to calculate the upper bound 
for Tc at Q = 0. Evaluating numerically the value of T such that J2j < 1, % > 0, 
we obtain, by sufficient condition (see Eq. (fT7|) ) , the upper bounds for as a function 
of Q, shown in figure [Tj together with the curve for the mean field results. We use 
fl23|) for the one-dimensional two-spin correlation function in obtaining figure [TJ This 
curve represents the rigorous upper bounds for T^, Qc- In particular, the mean field 
values are Qf^^ = zJ/2 and T*^^^(0) = zJ/ks- The values for and Tc{0), which 
are the rigorous upper bounds ioi d = 2, z = 4 and d = 3, z = 6, obtained in the 
present calculation are: (a) d=2, z=4; = 3.014 and Vic = 1.3755 J, (b) d=3, z=6; 
^ = 5.423 and Qc = 2.4466J. 

In table [T], we compare the results obtained by the effective field calculation (EFT) 
[0], the high temperature expansion (HTE) [3 E] and the present results for Qc- 




Figure 1. Phase diagram of the TIM model for z = 4 and z = 6. Solid line is the 
mean field approximation result. Results for our approach evaluated by using the spin 
correlation function separated by two sites, Eg. ([25)1 are given by the dotted {z = 4) 
and dashed {z — 6) lines. 
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. Estimativcs for 


n^/QcMFA for d 


=2 and d=3. 




d = 2,Z = 4: 


d = 3, Z = 6 


MFA 


1 


1 


EFT 


0.688 


0.784 


HTE 


0.770 


0.860 


Present work 


0.643 


0.813 



5. Concluding remarks 

In this paper we have obtained rigorous upper bounds for the critical couplings of the 
transverse Ising model. The procedure was based on an approximation for an exact 
identity for the two-spin correlation functions and on rigorous inequalities for the spin 
correlation functions. The approximated relation for the two-spin correlation function, 
Eq.(|5]), used in this procedure, is consistent with the rigorous inequalities, Eq. (fT7|l . since 
both act in the same direction of the inequahties. The upper bounds were applied for 
two- and three- dimensional models. 
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